Affine and Non-AfRne Motions in Sheared Polydisperse Jammed Emulsions 
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We study dense and highly polydisperse emulsions at droplet volume fractions (j) > 0.65. We apply 
oscillatory shear and observe droplet motion using confocal microscopy. The presence of droplets 
with sizes several times the mean size dramatically changes the motion of smaller droplets. Both 
afBne and nonafRne droplet motions are observed, with the more nonafBne motion exhibited by the 
smaller droplets which are pushed around by the larger droplets. Droplet motions are correlated 
over length scales from one to four times the mean droplet diameter, with larger length scales 
corresponding to higher strain amplitudes (up to strains of about 6%). 

PACS numbers: 83.80, 47.57.Qk, 83.85.Ei 
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Many materials behave as solids, despite having liquid- 
like amorphous structure. Obvious examples are glassy 
materials, which are disordered on an atomic scale 
Other examples include piles of sand, dense colloidal 
pastes, and shaving cream foams, which are disordered 
on the scale of microns or millimeters. Solids are char- 
acterized by their response to stress: they deform elas- 
tically (below a yield stress), rather than flowing. This 
behavior is well understood for crystalline materials, but 
somewhat more poorly for amorphous materials [l|, 
In particular, most work has focused on the case of fairly 
monodisperse samples composed of particles all of similar 
sizes [3l4l5|. However, many natural materials of inter- 
est are highly polydisperse, with particle sizes varying by 
factors of ten or more. The rheology of such materials 
has been much less widely studied [ll, • There have 
been differences noted from the monodisperse case, such 
as a lower strain amplitude required for viscous flow [Tsj 
and diminished sample viscosity [l^, H^l • The causes of 
these differences are not well understood. 

In this Letter, we study highly polydisperse emulsions. 
Our emulsions are composed of oil droplets in water, sta- 
bilized by a surfactant, and are at sufficiently high vol- 
ume fractions (d> > 0.65) that the samples act as amor- 
phous solids [13, [H, H^l ■ We subject the samples to low 
amplitude oscillatory shear and follow the droplet mo- 
tion in the interior of the sample via confocal microscopy 
[l3[lll,[2l|. Most droplets rearrange elastically Q and 
move sinusoidally. However, these motions are not nec- 
essarily affine, as shown in Fig. [T] In particular, our 
main finding is that in a highly polydisperse emulsion. 
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the smaller droplets frequently undergo highly nonaffine 
droplet motion. In contrast to the shear of monodis- 
perse samples [1, [13, IHl, , the largest droplets allow 
for "cross-talk" between layers at different heights which 
have different mean velocities. The motions of droplets 
are correlated over large length scales, up to four times 
the mean droplet radius, with longer range correlations 
found for higher applied strain amplitudes. 

We use the shear-rupture method of Ref. [13] to create 
decane-in-water emulsion droplets stabilized with SDS, 
skipping the fractionation step. The continuous phase is 
composed of a 65:35 volume ratio of water and glycerol 
to index match the decane droplets. Volume fractions 
are tuned to the range 0.65 < (/> < 0.85 by centrifuga- 
tion and dilution. Prior to our experiments, the samples 
are gently stirred to prevent any size segregation due to 
sedimentation, although for our high volume fractions 
sedimentation and size segregation are exceedingly slow. 

We place the samples in a parallel-plate shear cell. The 
design is similar to that of Ref. Q and our apparatus is 
described more fully in Ref. [ll|. The gap of the cell is 
fixed at H — 200 /im. The lower glass plate does not 
move, whereas the top plate is driven sinusoidally at a 
frequency / = 1 Hz by a piezoelectric driver. The ampli- 
tude is A = 40 /im and therefore the macroscopic strain 
amplitude is 7 = A/H = 0.2. Where the sample comes 
into contact with the glass plates, we add a coating of 
ScotchGard (3M). The sample wets the ScotchGard, pin- 
ning the sample to the coated region and ensuring a no- 
slip boundary condition at each plate. The oil droplets 
are deformable, but our experiments are conducted at 
volume fractions and strain amplitudes over which the 
droplets are fairly spherical. 

The sample is imaged from below through the station- 
ary plate using a confocal microscope. Fluorescein dye is 
added to the continuous phase to visualize the droplets as 
shown in Fig. [T] To quantify the size distribution of our 
system we acquire three-dimensional (3D) image stacks 



Figure 1: Confocal microscopy image of our polydisperse 
decane-in-water emulsion at (p — 0.65. The arrows indicate 
the total displacement of each droplet from the peak-to-peak 
of one half oscillation cycle; they are twice the length of the 
true displacement, for easier visualization. The mean droplet 
displacement (2.8 to the left during the half-cycle) has 
been subtracted off from the displacement arrows. The image 
corresponds in time to the mid-point of this half-cycle. The 
width of the image is 56 fim and the arrows are to scale. The 
depth is 2 = 24 fim, the strain amplitude is 7 = 0.067, and 
the volume fraction is rA = 0.65. 



from a static sample of size 56 x 59 x 80 /im"^ at a verti- 
cal resolution of dz = 0.1 fiui. To observe the dynamics 
when sheared, we must take data as rapidly as possible 
and thus we only use two-dimensional (2D) images. For 
the 2D experiments, images of size 56 x 59 /im^ are ac- 
quired at a rate of 90 images per second for 33 s. 

Using the 3D data from static samples, we determine 
the true droplet radii using custom software implement- 
ing the method of Ref. [231. The size distribution ob- 
tained from this method is shown in Fig. [51 The aver- 
age droplet radius is 1.2 /xm, the standard deviation is 
0.6 /Ltm, and the Sauter mean radius = {r^)/{r'^) is 
2.3 /Ltm. The data for small radii are reasonably well fit 
by an exponential with decay length 0.5 yum. However, 
there are a number of unusually large droplets with sizes 
larger than 5 fim. While these droplets are rare, they 
account for a nontrivial portion of the volume, as can 
be seen by the volume-weighted probability distribution 
shown in the inset to Fig. [21 

For 2D data analysis, we use a slightly different anal- 
ysis technique. We identify droplets using the 2D-Hough 
transform [24], which lets us identify the droplets' radii 
and positions in each image. From that data, we then 
use conventional techniques to track their corresponding 
trajectories [25|. Both standard tracking and the itera- 
tive image tracking technique described in Ref. [2l| are 
used to reconstruct each trajectory. Note that for each 



Figure 2: Probability distribution of droplet sizes in a = 
0.80 sample. Droplets with radius r < 5 ^m are fit to an 
exponential (dashed line) with decay length 0.5 nm. The data 
for r > 11 nm correspond to only three observed droplets, 
and thus the true shape of the distribution is ill-defined for 
the larger droplets. The inset shows the same data, with the 
probability weighted by the droplet volume by multiplying by 
r'^. Note that the inset is plotted on a linear scale, while the 
main graph is plotted on a logarithmic scale. 



droplet, because our observation is only in 2D, we do 
not know the true radius t^d. However, we observe that 
due to the curvature of the interfaces, droplets are not 
distinguishable in 2D slices when they are viewed more 
than Az w 0.6r3D away from their center, and so for each 
droplet, the droplet radius r we apparently observe is in 
the range O.SrsD < r < rso. 

When an oscillatory strain is applied to the emulsion, 
the majority of droplets rearrange reversibly and peri- 
odically at the driving frequency. The average velocity 
field is v{z,t) — 72: sin(27r/i)a;. The largest droplets 
^ 8 ^m) move sinusoidally at this mean velocity; 
this is reflected by the very short displacement vectors 
drawn on the largest droplets in Fig. [1] In contrast, 
smaller droplets move in a variety of directions. Some 
move affinely with the applied shear, but many move 
in other directions, as shown in Fig. [H This behavior 
is likely due to the largest droplets. At the equator of 
a large droplet, it moves with the expected velocity for 
that height z. However, droplets tend to move coherently 
as spheres (rather than deforming) , and thus the top of a 
large droplet moves with an amplitude jr too small rel- 
ative to the expected velocity at height z + r. Likewise, 
at height z — r the large droplet is moving more quickly 
than the mean velocity for that height. Thus, two large 
droplets that are nearby but with centers at different z 
values do not have the same velocities, and as they move 
back and forth sinusoidally, they push and pull on the 
smaller droplets between them. These smaller droplets 
thus have apparently random motions. In practice, the 
larger droplets are rarer and thus less likely to influence 
each other. Moreover, they move based on the average 
influence of the smaller droplets surrounding them, and 
so their motion is likely to tend to average velocity fleld. 
The contrast in motion between large and small droplets 
in a polydisperse sample is the main conclusion of this 
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Table I: Volume fraction 0, observed local strain amplitudes 
7ioc for different volume fractions (j>; observed slip length zq; 
and the characteristic length scales and ^y. The macro- 
scopic applied strain amplitude is 0.20 for all experiments. 
The uncertainties for all listed lengths are ±0.3 /im. 



Letter. This differs qualitatively from cases where large- 
scale flows cause non-afRne motion, and which typical 
require large amplitude strain 0, [l3| . 

Over the 33 s movies, approximately 8% of the droplets 
make an irreversible rearrangement at some point. This 
only occurs with smaller droplets, r <^ 5 /xm. Before 
and after the irreversible rearrangement, the droplets 
move periodically. The rarity of plastic rearrangements 
in our data is similar to a prior study of a more monodis- 
perse emulsion [sf. 

At this point we only study droplets whose trajectories 
are reversible (elastic) and thus periodic. We respectively 
denote as x{t) and y{t) the components of a trajectory 
parallel and perpendicular to the shear axis at time t. 
A least squares fit is applied to each component with 
functional forms: 
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Figure 3: Probability distributions of the fit parameters. 
These are (a) amplitudes of the parallel components, (b) am- 
plitudes of the perpendicular components, (c) phase angles 
of the parallel components and (d) phase angles of the per- 
pendicular components of the trajectories. Data correspond 
to a, cj> = 0.65 emulsion at a local oscillatory strain amplitude 
7 = 0.070, at a depth z = 2A /im. The data for P(ax) are cen- 
tered around {a^) = 2.4 /im. The means for the phase angles 
are arbitrary as they depend on when we set t = 0, although 
note that the distribution for 9y is centered around the mean 
value of Ox- The lack of symmetry in these distributions is 
due to observing a finite number of droplets. 



x{t) = ax sin(cL't + Ox) 

y(t) = aySm{u!t + 9y) (1) 

where cj = 27r/, using the known value of the driving fre- 
quency /. The fit parameters ax and ay give the ampli- 
tudes of the periodic motions of the droplet, and 9x and 
9y are phases. These functional forms provide a good fit 
to the particle trajectories. 

We study the sample at depths z — 24,36,48, and 
60 /im, relative to z = defined at the stationary bottom 
plate. At each height we compute the mean amplitude 
(ax)- We find (ax) ^ 7ioc ■ [z — zq) where 710c is the local 
strain amplitude and zq is a slip length, (ax) does not 
extrapolate to at z = 0, but rather at negative values 
ranging from zq = —8 to —14 /im. The local strain is al- 
ways smaller than the applied strain. This suggests that 
the emulsion partially slips at the top plate, or possibly 
has a shear band somewhere where the local strain is 
significantly higher. See Table I for further information 
about each particular experiment. We do not observe 
the behavior of purely slipping emulsions in any of our 
experiments [2d| . 

The distributions of the fitting parameters ax, ay, Ox 
and Oy are quite broad, as shown in Fig. [31 While many 
droplets move with the mean amplitude {ax) as appro- 
priate for that height, several have amplitudes that dif- 
fer by 0.5 /xm or more from the mean. Negative values 
of ax — {ax) indicate droplets moving with smaller am- 
plitudes than might be expected, and likewise positive 
values indicate droplets moving with larger amplitudes. 



These results are equally true in the direction of the ap- 
plied strain [Fig. Efa)] and perpendicular to this direc- 
tion [Fig. El^b)], showing many droplets have significant 
nonaffine motion. The amplitudes associated with big 
droplets are found at the central peaks of P{ax) and 
P{ay). The outliers are more likely to be associated 
with smaller sized droplets, with the tails of the distri- 
bution coming from some of the smallest droplets with 
r ^ 1 /im. Note that ay from our fits (Eqns.[T|) is positive: 
to get values Oy < 0, we assume that all droplets move in 
phase, and so droplets with phase angles Oy that appear 
TT out of phase with the dominant motion are adjusted, 
ay —o-y, Oy — 7> {Oy — tt). In general we find {ay) « 0, as 
expected by symmetry. 

The phase angles also vary from droplet to droplet. 
In the direction parallel to the shear, most droplets 
move in phase, as P{9x) is peaked and narrow as shown 
in Fig. ^c). In contrast, P{Oy) is broadly distributed 
around a central peak, as shown in Fig. [3]Jd). The 
phase angle in y ranges up to 7r/2 away from the av- 
erage phase, indicating some droplets move periodically 
but completely out of phase with the driving strain. 

To understand the spatial character of the particle be- 
havior, we construct images with each droplet positioned 
at its time-averaged position, and color coded based on 
its motion, shown in Fig. [H The color scale corresponds 
to the value of a^, (left) and ay (right). The biggest 
droplet has a color corresponding to the center of the 
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Figure 4: (Color online) Droplets which move elastically, 
drawn at their mean position. The color code indicates the 
parallel (left picture) and perpendicular (right picture) am- 
plitude of each droplet. The color bar denotes the amplitude 
scale in /im, where denotes the average amplitude for these 
data {{ax) = 2.4 fj,m, and using for the ay data). The 
images correspond to the data shown in Fig. [S] 



color scale, in other words, it is moving with the aver- 
age motion of all droplets. The droplets moving with 
anomalous amplitudes corresponding to the extremes of 
the color scale are the smaller droplets. Droplets with 
similar Ux or ay tend to be close together. 

To gain further insight into the spatial clustering of 
droplets with similar characteristics, we identify droplet 
pairs which are separated by a surface-to-surface distance 
Ar, that is, droplets separated by a center-to-center dis- 
tance ri2 = ri+r2 + Ar, where is the radius of droplet 
i in the 2D image. Then, we compute the Pearson corre- 
lation coefficient of the amplitude ax '■ 



C{Ar) 



1 



{ax, - {ax)){axj - {ax}) 



(2) 



where N{Ar) is the number of neighboring droplets, and 
axi and a^^ are the x-amplitudes of droplets i and j. 
The average (ax) is for all droplets comprised by N{Ar) 
and cTa^ corresponds to the variance of the distribution 
of ax- The choice of using the surface-to-surface dis- 
tances rather than center-to-center is perhaps not ob- 
vious. However, each individual droplet moves as a solid 
object, thus completely correlated with itself (distances 
1^12 < ''i)- Examining the surface-to-surface motion lets 
us avoid considering this artificially correlated solid-body 
motion, and instead probe the properties of the effective 
medium between droplets. If we consider center-to-center 
separations, the results are noisy and do not have a sim- 
ple dependence on the distance. 

Figure EJa) shows these correlation functions for ax 
(solid symbols) and ay (open symbols), for one sample 
at two different strain amplitudes. The correlation func- 
tions exhibit exponential decay with decay lengths in the 
range of 8-15 /xm; the decay lengths for different samples 
are listed in Table L These lengths are comparable to 
the sizes of the larger droplets in the sample. Figure[5lja) 
shows that correlations in the larger strain case (squares) 
decay slower than than in the small strain case (circles). 



We find £,x and ^j, depend on 710c as shown in Fig.[5fb) for 
this (f) — 0.65 sample. Given the dependence of £ on 7iocj 
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Figure 5: (Color online) (a) Spatial correlation functions 
of the droplet amplitudes a^ (solid symbols) and ay (open 
symbols). The strain amplitudes are 710c = 0.017 (circles) 
and 710C = 0.085 (squares). The correlation functions are 
normalized by their value at Aro = 2 /im. The data are 
from the same sample as Figs. [3] and [H (b) Values of the 
exponential decay length scale for this sample as a function 
of strain amplitude, for correlations of a^ (solid symbols) and 
ay (open symbols). The uncertainty of 7 is 6%. The error bars 
are the standard deviations found from multiple experiments. 
All data for both panels of this figure are taken at depth 
= 24 /im. The lines are guides to the eye. 



it suggests that perhaps the variation we see for samples 
with different is perhaps more due to the variability in 
710C for those data; see Table I. 

We have studied jammed polydisperse emulsions, and 
observed highly complex droplet motion when our sam- 
ples are sinusoidally sheared. Most droplets move peri- 
odically, but different droplets have different amplitudes 
and phases of their motion. The large droplets greatly 
rearrange the small droplets, although nearly all of this 
motion is reversible. We find length scales over which 
the droplet motions are correlated. These length scales 
range from 1 to 4 times the average droplet diameter, 
with the largest values found for the highest strains (up 
to strains of ~ 0.08). Overall, our results suggest that 
the rheological response of highly polydisperse systems 
has a complex microscopic origin: larger droplets cause 
long-range correlations in the dynamics, while smaller 
droplets are like an effective fluid. Theoretical descrip- 
tions derived for less polydisperse systems will likely not 
apply or need to be modified [1, @] ■ 
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